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Equilibrium phase transitions may be denned as nonanalytic points of thermodynamic functions, 
e. g., of the canonical free energy. Given a certain physical system, it is of interest to understand 
which properties of the system account for the presence, or the absence, of a phase transition, 
and an investigation of these properties may lead to a deeper understanding of the physical 
phenomenon. One possible way to approach this problem, reviewed and discussed in the present 
paper, is the study of topology changes in configuration space which, remarkably, are found to 
be related to equilibrium phase transitions in classical statistical mechanical systems. For the 
study of configuration space topology, one considers the subsets M v , consisting of all points 
from configuration space with a potential energy per particle equal to or less than a given v. 
For finite systems, topology changes of M v are intimately related to nonanalytic points of the 
microcanonical entropy (which, as a surprise to many, do exist). In the thermodynamic limit, 
a more complex relation between nonanalytic points of thermodynamic functions (i.e., phase 
transitions) and topology changes is observed. For some class of short-range systems, a topology 
change of the M v at v = vt was proven to be necessary, but not sufficient, for a phase transition 
to take place at a potential energy v%. In contrast, phase transitions in systems with long- 
range interactions or in systems with nonconfining potentials need not be accompanied by such 
a topology change. Instead, for such systems the nonanalytic point in a thermodynamic function 
is found to have some maximization procedure at its origin. These results may foster insight into 
the mechanisms which lead to the occurrence of a phase transition, and thus may help to explore 
the origin of this physical phenomenon. 
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Preface 

It was at the end of the 1990s when the application of 
concepts from differential geometry to Hamiltonian dy- 
namical systems led to a conjectured connection between 
the occurrence of equilibrium phase transitions in classi- 
cal Hamiltonian systems and some topological quantities 
of configuration space subsets of these systems. Since 
then, the interest in this approach and the number of 
people working on the topic has increased, and so has 
the number of results. 

At the time of this writing, an overview of the sub- 
ject is difficult to attain. First, the results are scattered 
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among a considerable number of publications and, sec- 
ond, several of the results, although correct, demand a 
reinterprctation as a consequence of recent findings and 
developments. The purpose of the present paper is to 
assemble from the known results, as far as possible, a co- 
herent picture of the relation between phase transitions 
and configuration space topology, and to indicate new 
lines of research which might open up from these con- 
cepts. 



I. INTRODUCTION 

Phase transitions, like the boiling and evaporating of 
water at a certain temperature and pressure, are com- 
mon phenomena both in everyday life and in almost any 
branch of physics. Loosely speaking, a phase transition 
brings about a sudden change of the macroscopic proper- 
ties of a many-particle system while smoothly varying a 
parameter (the temperature or the pressure in the above 
example). Probably the main reason for the unabated 
interest that phase transitions have received already for 
more than a century is their omnipresence in all branches 
of physics (and also in related fields like biology or engi- 
neering): be it the formation of stars in astrophysics, the 
transition to superconductivity in solid state physics, or 
the opening of the DNA helix in biophysics, examples of 
many-particle systems which undergo a phase transition 
are widespread and of indisputable relevance in science. 

Phase transitions can occur in both equilibrium and 
noncquilibrium systems, but the focus will be exclusively 
on the equilibrium case in this exposition. There is a 
plethora of books on the subject, ranging from experi- 
mental to theoretical and mathematical treatises. Espe- 
cially the theory books are to a large extent concerned 
with the classification of different types of phase tran- 
sition: remarkably, very different physical systems may 
show quantitatively identical properties in the vicinity of 
the phase transition point, and this fascinating topic of 
universality, explained using the renormalization group 
theory, has attra cted a lot of interes t espec i ally in the 
1970s and 1980s (jBinnev et all Il992t iLesnd . Il998l ). In 
the present paper, instead of discussing the characteris- 
tics of the different types of phase transition, we take one 
step back and inquire about the conditions which may 
lead to the occurrence of a phase transition in a given 
system. 

The mathematical description of equilibrium phase 
transitions is conventionally based either on Gibbs mea- 
sures on phase space or on (grand)canonical thermody- 
namic functions, relating their loss of analyticity (or, in 
other words, the appearance of a singularity) to the oc- 
currence of a phase transition. 1 A nonanalyticity of a 



thermodynamic function leads to a discontinuity or to a 
divergence in some derivative of this function, and this 
is a typical signature of a phase transition as measured 
experimentally. Within the (grand)canonical ensemble of 
statistical mechanics, such a nonanalytic behavior can oc- 
cur only in the thermodynamic limit, in which the num- 
ber of degrees of freedom TV of the system goes to infin- 
ity. 2 Conceptually, the necessity of the thermodynamic 
limit is an objectionable feature: first, the number of de- 
grees of freedom in real systems, although possibly large, 
is finite, and, second, for systems with long-range inter- 
actions, the thermodynamic limit may even be not well 
defined. These observations indicate that the theoretical 
description of phase transitions, although very successful 
in certain aspects, may not be completely satisfactory. 

Apart from this conceptual shortcoming, in the field 
of phase transitions there arc many problems of applied 
nature which are far from being settled. One of those 
is the search for sufficient or necessary conditions for 
the occurrence of a phase transition. Among the neces- 
sary conditions for the occurrence of a phase transition, 
there are some of reasonable generality, lik e the Mermin 



Wagner theorem and its general i zations dFannes et al 
1984 iFrohlich and Pfisterl . 1 198 ll iMermin and Wagner 



19661 ) or the theorems on the absence of p hase tran- 
sitions in cert ain one-dimensional system s by Ivan Hove 
(195(3) and bv lCuesta and Sanchez! ( 2004 ). Yet improved 
criteria are of course desirable. Less is known about 
conditions sufficient to guarantee a p hase transition to 
take place. The Peierls argume nt dPeierlsl. 19361) or 
the Frohlich-Simon-Spencer bound (jFrohlich et 
can be used to prove the existence of phase transitions 
without explicitly computing a thermodynamic poten- 
tial, but their application is model specific and may be 
difficult depending on the system of interest. 

The above considerations motivate a further study 
on the "nature" of phase transitions, of the underlying 
mechanisms leading to a nonanalytic point of a thermo- 
dynamic function, and of the conditions under which they 
can occur. A classic result identifying such a nonana- 
lyticity generating mechanism is the seminal theorem of 
Lee and Yang, relating the properties of the zeros of the 
grandcanonical partition function in the complex fugacity 
plane to nonana lyticities of the corre sponding thermody- 
namic function ( Lee and Yangl . Il952) . The main issue of 
the present paper is to investigate the mechanism which 
is at the basis of a phase transition using a different ap- 
proach, based on concepts from differential geometry and 
topology. 

This topological approach emerged from the study 
of Hamiltonian dynamical systems and is therefore — at 
least in its present formulation — applicable to classical 
(i.e., nonquantum mechanical) systems. Hamiltonian 



1 Standard refer ences in math emat ical physics for these two points 
of view are bv lRuelld <196Sl) and lGeorgiil Jl98Sf) . 



2 Apparently it was first suggested by H. A. Kramers in 1937 to 
take this limit. 



dynamics can be viewed as a geodesic flow on the configu- 
ration space, provided the latter is equipped with a suit- 
able metric. 3 By numerical methods, geometric quan- 
tities of such metric spaces of Hamiltonian dynamical 
systems were studied as a functio n of energy o r tem- 
pera ture in a serie s of papers [see ICasetti et al\ ( 2000f ) 
and iPettinil (|2007l ) for reviews]. For a system under- 
going a phase transition in the thermodynamic limit, 
these geometric quantities display discontinuous or cusp- 
like features remarkably close to the transition energy 
or temperature. One possible mechanism behind such 
a dramatic change of geometric quantities of configura- 
tion space can be a change of its topology, and it was 
this line of reasoning which led Pettini and co-workers 
to conjecture a connection between phase transitions, on 
the one hand, and topology changes within a family of 
certain configuratio n space subsets, on the other hand 
( Caiani et al\ . Il997t ). Subsequently, such a connection 



was proven to hold true for a certain class of systems with 
short-range interactions, showing that a topology change 
is a n ecessary condition fo r a phase transition to take 
place (jFranzosi and Pettinil . |2004 l2007t iFranzosi et all 
120071 ). This theorem, for the class of systems covered 
by its assumptions, suggests the interpretation of certain 
topology changes as the relevant mechanism behind the 
generation of a nonanalyticity in a thermodynamic func- 
tion, and it furthermore allows one to exclude the occur- 
rence of a phase transition when such topology changes 
are absent. 

The use of concepts from topology to describe a physi- 
cal phenomenon is particularly appealing due to the fact 
that topology yields a very reductional description: con- 
sidering only the topology of, say, a surface, a significant 
amount of "information" (on curvatures, for example) 
is disregarded, and only a small part (like connectivity 
properties) is kept. If one then succeeds to capture the 
essentials of the phenomenon of interest with the remain- 
ing information only, the resulting description will be an 
efficient one, and one might hope to get an unblurrcd 
view onto the mechanism which is at the basis of the 
phenomenon. 

The initial hope that such a topological approach 
might provide a unified and completely general descrip- 
tion of phase transitions turned out to be over-optimistic. 
It was only in the last few years that evidence accumu- 
lated, disproving the general validity of the hypothesized 
connection betwe en phase transi t ions and configuration 
space topology ( Angelani et all l2005bt iBaroni 



3 Both Hamiltonian dynamics and geodesies of a Riemannian man- 
ifold can be denned by some variational principle: the trajecto- 
ries of Hamiltonian dynamics are the extrema of the Hamiltonian 
action functional, whereas the geodesies of any Riemannian man- 
ifold are given by the extrema of the arclength functional. It is 
this structural similarity which al lows for a geometric formu la- 
tion of Hamil t onian dynamics [see iMarsden and Ratiu ( 1 9941 ) or 
ICasetti et oil l l2000ri l. 



iGaranin et al. l, l2004HHahn and Kastneil 12005). This ob- 
servation, perceived as a major set-back at first, ap- 
pears less dramatic in light of subsequent findings. How- 
ever, and this was one of the motivations for writing the 
present paper, it alters the understanding of the topolog- 
ical approach to phase transitions, as well as of several 
of the results of model calculations reported in the liter- 
ature. The aim of the present paper is to review results 
on the relation between configuration space topology and 
analyticity properties of thermodynamic functions. This 
relation is found to depend on the physical situation, in 
particular on whether the system of interest is finite or 
infinite and whether the interparticle interactions are of 
short range or of long range. The results may help to 
deepen the understanding of the basic mechanisms be- 
hind phase transitions in the infinite system case, and 
they can explain the peculiar (non)analyticity properties 
of microcanonical entropy functions of finite systems. 

The paper is structured as follows: We start by fixing 
notations and giving some basic definitions used in the 
topological approach to phase transitions in Sec.|TIJ Since 
Morse theory provides a suitable mathematical frame- 
work for the study of the topology of the configuration 
space subsets of interest, a summary of elementary re- 
sults of this theory as well as an application to a sta- 
tistical mechanical model is given in Sec. IIII1 The re- 
lation between certain topology changes and nonanalyt- 
icities in thermodynamic functions of finite systems is 
discussed in Sec. IIVI The rest of the paper is devoted 
to the more intricate case of infinite systems. A theo- 
rem rigorously establishing a connection between phase 
transitions and topology changes in configuration space 
is presented in Sec. LY1 where also results on the config- 
uration space topology of models beyond this theorem's 
assumptions are reviewed. In Sec. I VII the limitations of 
the topological approach are explored by studying two 
models for which the proposed relation between phase 
transitions and configuration space topology turns out 
to be invalid. Some proposed sufficiency conditions on 
the topology changes, guaranteeing the occurrence of a 
phase transition, are critically discussed in Sec. I VIII We 
conclude the paper with a summary in Sec. IVIIII 



II. DEFINITIONS AND PRELIMINARIES 

Discussing topology and topology changes in the Intro- 
duction, vague reference was made to certain subsets of 
configuration space. It is the aim of the present section 
to fix notations and to define configuration space subsets 
and thermodynamic quantities to which reference will be 
made throughout. 



A. Standard Hamiltonian systems 

We consider classical Hamiltonian systems consisting 
of N degrees of freedom, characterized by some Hamil- 
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tonian function 

H : A^v — > H (2.1) 

which maps the phase space A at C R 2JV onto the reals. 
For convenience we assume H to be of the standard form 

H(p;q) = ±(p,p) + V(q), (2.2) 

where p = (pi, ■ . . , jJjv) i s the vector of momenta and 
g = (qi, • • • ,Qn) is the vector of position coordinates. 
Masses have been set to unity and (•, •) denotes the usual 
scalar product. The potential 

^:rV->lR (2.3) 

maps the configuration space rV C WL N onto the reals. 
The restriction to Hamiltonian functions of the standard 
form (12. 2p is convenient, but not essential for the con- 
cepts discussed. Throughout the paper, we assume rV 
to be continuous (in contrast to classical spin models like 
the Ising model or the Potts model which have a discrete 
configuration space) . This will be important for the kind 
of configuration space topology we consider in the follow- 
ing. 



B. Configuration space subsets 

From the potential function V, we define the family of 
subsets {■M„}, ug]ft , where 

M v = V~ l (-co, vN] = {qtT N \ V(q) < vN} . (2.4) 

V" 1 gives the preimage of a set under V, hence A4 V is 
the subset of all points q from configuration space Tjv for 
which the potential energy per degree of freedom V(q) /N 
is equal to or less than a given value v. Similarly, the 
related family {S„} lJgR can be defined, where 

= V~ l (vN) = {qeT N \ V(q) = vN} (2.5) 

consists of all points q from the configuration space 
for which the potential energy per degree of freedom 



4 For the thermodynamics of phase transitions, a quadratic form 
in the momenta as i n Eq. I|2.2|| only leads to a shift in the free 
energy [however, see lCasetti and Kastnen l|2007f ) for a discussion 
of the pitfalls of this reasoning]. Furthermore, after becoming 
familiar with the topological concepts and notations, the reader 
may convince himself that any critical point q c of the potential 
V corresponds to a critical point (0; q c ) of a standard Hamilto- 
nian function, and that therefore the contribution of a standard 
kinetic energy as in Eq. (12.21) to the topological approach is a 
trivial one. As a consequence, the results presented remain valid 
also for models without a kinetic energy term (which is the situ- 
ation typically encountered when studying spin systems like the 
Heisenberg model). 



V(q)/N is equal to a given value v. These constant po- 
tential energy subsets form the boundaries of the corre- 
sponding M. v , i. e., 

E« = dM v , (2.6) 

so that A4 V or £„ are closely related. It is therefore a 
matter of convenience to use one quantity or the other, 
depending on the actual situation of interest. The main 
topic of this paper is the relation of topology changes 
of A4 V or Y. v to nonanalytic points of thermodynamic 
functions. 



C. Thermodynamic functions 

Strictly speaking, the notion of a thermodynamic func- 
tion should be restricted to functions describing the equi- 
librium behavior of systems in the thermodynamic limit 
of infinitely many degrees of freedom. For matters of con- 
venience, we likewise will speak of thermodynamic func- 
tions when referring to their finite system counterparts. 

Regarding their analyticity properties, thermodynamic 
functions obtained from different statistical ensembles 
can differ drastically. In this paper we discuss analyt- 
icity properties of microcanonical and canonical thermo- 
dynamic functions. 5 

1. Microcanonical thermodynamic functions 

The microcanonical ensemble provides the framework 
for the statistical description of an isolated physical sys- 
tem in which the total energy is conserved. The funda- 
mental quantity of this ensemble is the Boltzmann en- 
tropy or microcanonical entropy as a function of the en- 
ergy (per degree of freedom) e, 

s JV ( £ ) = lln/ dpdq6[H(p;q)-Ne], (2.7) 
iV Ja n 

where 6 denotes the Dirac distribution. 6 A related quan- 
tity is the configurational microcanonical entropy as a 
function of the potential energy (per degree of freedom) 
v, 

s N (v) = 4 In / dqS [V(q) - Nv] . (2.8) 

2. Canonical thermodynamic functions 

The canonical ensemble provides the framework for the 
statistical description of a system coupled to an infinitely 



5 The extension to further statistical ensembles is straightforward. 

6 Here and in the following we define thermodynamic functions 
always per degree of freedom, which accounts for the factor 1/JV 
in the definitions. 
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large heat bath of inverse temperature (3. The fundamen- 
tal quantity of this ensemble is the canonical free energy 

= — ik^f dpdge-" ff fe«). (2.9) 
Np J An 

A related quantity is the configurational canonical free 
energy 

SM= ~w ln L dqe ~ mq) - (2 - 10) 



3. Relation of microcanonical and canonical thermodynamic 
functions 



It follows from the definitions l|2.7j) and (|2.9|) that mi- 
crocanonical entropy and canonical free energy are re- 
lated by 

f N (® = -lk ln I A ^ N£ ) ew{N[s N (e)-pe}}. 

(2.11) 

Similarly, for the configurational thermodynamic func- 
tions the relation 

M3) = -j^k / d ( Nv ) exp {N [s N (v) — /3v]} 

(2.12) 

holds. In the thermodynamic limit, using Laplace's 
method for the evaluat ion of asymptotic integrals 
([Bender and Orszad. fl999h , well-known relations of ther- 
modynamic functions by means of Legendre-Fenchel 
transformations, 



Pfoo(P) = sup [soo(e) - (3e] 



(2.13) 



and 



-MM = sup [««,(«) (2.14) 
arc obtained from Eqs. (|2~TT|) and (j2T2ll . 

D. Nonanalytic points and phase transitions 

The main topic of the present paper is the relation 
between topology changes within the families {M v } v€R 
or {Ti v } v& - SR of configuration space subsets defined in 
Sec. III.B1 on the one hand, and nonanalytic points of 
thermodynamic functions, on the other hand. 

Definition II. 1. A nonanalytic points is a point in the 
interior of the domain of a real function at which the 
function is not infinitely many times real-diffcrentiablc. 
Synonymously we speak of nonanalyticities of the func- 
tion. 

Different definitions of thermodynamic phase transi- 
tions can be found in the literature, where the most 



common ones are based cither on the (non)uniqucness of 
translationally invariant Gibbs measures on phase space 
or on the (non)analyticity of thermodynamic functions. 7 
Throughout this paper we will use the following version 
of the latter approach. 

Definition II. 2. A phase transition is defined as a non- 
analytic point of the canonical free energy /jy. The tran- 
sition is called discontinuous if the first derivative of /jy 
is discontinuous, otherwise it is called continuous. 

In an older but somewhat misleading terminology, dis- 
continuous phase transition are called first-order phase 
transitions whereas continuous ones are referred to as 
second-order phase transitions. 

Computing first and second derivatives of the canoni- 
cal free energy, one arrives at experimentally measurable 
quantities like the canonical caloric curve 



u(f3) = -^/3f N ([3) 



or the canonical specific heat 



d_ 

d[3 



d(3 



(2.15) 



(2.16) 



and the effect of the "order" of the phase transition is 
enhanced in the behavior of such derivatives. For exam- 
ple, in the presence of a continuous phase transition, the 
caloric curve u(/3) shows a kink at the transition tem- 
perature p\, whereas a discontinuity in u is observed in 
the case of a discontinuous phase transition (hence the 
name). The liquid-vapor transition of water is an ex- 
ample of a discontinuous phase transition, whereas the 
transition from ferromagnetic to paramagnetic behavior 
is a continuous one for many materials. 

Note that definition III. 21 of a phase transition and its 
order explicitly refers to a particular statistical ensemble, 
namely the canonical one. This situation seems some- 
what unsatisfactory and has been the cause of a long- 
standing debate. Since, as mentioned in Sec. III. CI the 
analyticity properties of thermodynamic functions may 
depend drastically on the statistical ensembles used, the 
consequences of this choice of definition are considerable. 
One can show that the canonical free energy fjq of a sys- 
tem of N degrees of freedom is a smooth function for all 
finite N, and hence no phase transitions occur in finite 
systems (lGriffithstll972h . This is in contrast to the be- 
havior of microcanonical thermodynamic functions like 
the entropy which can have nonanalytic points also in 
the case of a finite number N of degrees of freedom. We 
discuss this issue in Sec. IIVI 



7 For some examples like the Ising model both definitions are 
known to coincide, for many others coincidence may be expected. 
However, counterexamples can be found as well, like the case of 
the Kosterlitz-Thouless phase transition of the two-dimensional 
XY model which has a u nique translationally i nvariant Gibbs 
state for all temperatures dBricmo nt et ni l Il977f). An introduc- 
tory discussion has been given by lLebowitj ||1999T ). 
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III. COMPUTATION OF TOPOLOGICAL QUANTITIES 

Computations of quantities characterizing the topol- 
ogy of the subsets M. v as defined in Eq. ()2.4|) can be 
found in the literature for several statistical physics mod- 
els. 8 A nume r ical c omputation has been reported by 
iFranzosi et al. (2000), whereas all other results are an- 
alytic, using methods from Morse theory in most of the 
cases. Morse theory plays an important role not only for 
the study of model systems, but also for general investi- 
gations on the relation between topology and nonanalyt- 
icities of thermodynamic functions. For these reasons, a 
short summary of some basic facts and concepts of this 
theory is given. 



Theorem III. 3. If g has no critical values in the interval 
[a, 6], then M a and Mb are homcomorphic, i.e., there 
exists a homeomorphism 10 mapping M a onto Mb- 

Homcomorphicity is synonymous to topological equiva- 
lence, so M a and Mb are topologically equivalent, M a ~ 
Mb, under the above stated conditions: no topology 
changes take place within the family {Mt} te r a b \ upon 
variation of the parameter t in the interval [a, b]. 



Proof of Theorem [lII.3[ see Chap. 
(120021) . 



3.1 of IMatsumotol 
□ 



At least in its standard form, Morse theory applies to 
the class of so-called Morse functions. 



A. Morse theory 

Morse theory establishes a link between the two mathe- 
matical disciplines of topology and analysis. For an intro- 
ductio n to Morse theory, see the tex t book by lMatsumotol 
(12002ft or the classic text bv lMilnorl (|l963f ). For the type 
of problem we are interested in, Morse theory allows to 
characterize the topology of the configuration space sub- 
sets Ai v defined in Eq. (|2.4j) by studying, by means of 
elementary analysis, the critical points 9 of the potential 
V. 

We consider a smooth (i. e., infinitely many times dif- 
ferentiable) function g : M — > M, mapping some Tri- 
dimensional manifold M onto the set of real numbers. 
In the context of our exposition, the role of this general 
function g will be played by the potential V of a standard 
Hamiltonian system as defined in Sec. III. Al 

Definition III.l. A point q c <G Al is called a critical 
point of g if the differential dg (q c ) at q c G M vanishes. 

Definition III. 2. A real number g c is a critical value of 
g if g (q c ) = g c for some critical point q c of g. 

With these definitions, we can state a first theorem 
relating properties of the function g to the topology of 
the subsets 



M t = {?£ M | g(q) t} 



(3.1) 



[defined analogous to the configuration space subsets A4 V 
in Eq. (l2~4l) ]. 



Such mo del calculation s were re ported by [Andronico et ql\ 
[12004); lAngelani et oil d2003l. l2005aT); iBaroml (|2002h 
Casetti et alj__ Jl999l. 120021 . |2003); IFranzosi et all [boom 

iGrinza and Mossal (|2004h 
Kast ner and Schnetj (12 006); 
1120041 ): Risau-Gusman et ali 



Gara nin et al\ 
Kastnerl [12004 



Ribciro Tcixcir a and Stariold 
1120051) . 

9 Not to be confused with critical points in the theory of phase 
transitions and critical phenomena. 



Definition III. 4. A critical point q c of g is nondegener- 
ate if the determinant of the Hessian of g at q c is nonzero. 

Definition III. 5. A function g : M — ► ]R is called a 
Morse function if every critical point of g is nondegener- 
ate. 

Then Morse theory relates the topology of M t to the 
critical points q c of g and their indices. 

Definition III. 6. The index of a critical point q c of g 
is the number of negative eigenvalues of the Hessian of g 
at q c . 

Theorem III. 7. If the interval [a, 6] contains a single 
critical value of g with a single critical point q c , then 
the topology of Mb differs from the topology of M a in 
a way which is determined by the index i of the critical 
point: Mb is homeomorphic to the manifold obtained 
from attaching to M a an i-handle, i. e., the direct product 
of an i-disk and an (to — z)-disk. 

Proof. A proof of this theorem, together with a pre- 
cise definit ion of "atta c hing a handle," can be found in 
Sec. 3.1 of IMatsumotol (|2002t) . A generalization to crit- 
ical values with more than one critical point is straight- 
forward and involves the attachment of more than one 
handle. □ 



With Theorems 1111 . 31 and IIII.7I we have transformed 
the problem of determining the topology of Mt to the 
problem of determining the critical points and critical 
indices of the underlying function g, which brings us back 
from topology onto the familiar grounds of analysis. 

At least in this standard version, the results of Morse 
theory apply only to the class of Morse functions speci- 
fied in Definition IIII. 51 Conceptually, this is an insignif- 
icant restriction, since Morse functions on M form an 
ope n dense subset o f the space of smooth functions on 
M (|Demazurd . [20001) . This means that, if the potential 



10 A homeomorphism is a continuous bijection between manifolds 
with continuous inverse. 
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V of the Hamiltonian system we are interested in is not 
a Morse function, we can transform it into a Morse func- 
tion V by adding an arbitrarily small perturbation, for 
example, 



N 



V(q) = V(q) + J2 hiQi 



(3.2) 



with some small hi £ R (i = 1, . . . , N). For practical 
purposes, however, adding a perturbation — and thereby 
destroying a symmetry present in V — may render the 
explicit computation of critical points and indices much 
more complicated (or even impossible). 

Theorem IIII.7I asserts that the topology of the man- 
ifolds M t is characterized by the critical points and 
their indices, but the "information" contained in these 
quantities — especially in the case of high dimensional 
manifolds — may be somewhat difficult to handle. A sim- 
pler, nonetheless useful characterization of the topology 
is given by the Euler characteristic, which can be ex- 
pressed by means of Morse numbers [see any textbook 
on algebraic topology for a definition of the Euler charac- 
terist i c for quite general topological spaces, for example, 
IVickl (119941) ]. 



Definition III. 8. The Morse numbers fa (i = 1,. .., m) 

of a function g on an m-dimcnsional manifold (or m- 
manifold) M are defined as the numbers of critical points 
of g with index i. 

Theorem III. 9. The Euler characteristic x of an m- 
manifold M can be expressed as the alternating sum of 
the Morse numbers fa of any Morse function on M, 



(3.3) 



i=0 



Proof. See, for example, Chap. 2, Sec. 9.3 of iFomenkd 
(119871) . □ 



An important property of the Euler characteristic x is 
that it is a topological invariant, i.e., different values of 
X for manifolds Mi and imply that M\ and are 
topologically nonequivalent. Hence monitoring the Eu- 
ler characteristic of the family {.M,,},,^ of configuration 
space subsets of some Hamiltonian system under varia- 
tion of the parameter u, we may get an impression of the 
way the topology of the A4 V changes. 



B. Model calculation: Mean-field fc-trigonometric model 

In an explicit calculation, depending on the methods 
applied (or applicable), different quantities characteriz- 
ing the topology of the subsets M v may be obtained. 
Typical examples are as follows: 

(i) the critical points q c of V and their indices (Defi- 
nitions HTT7T1 and |HT6|) ; 



(ii) the critical values v c — V(q c )/N at which the topol- 
ogy of the M. v changes (Definition IIII.2|) : and 



(iii) the Morse numbers fa of the M v (Definition IIII.8|I . 
which, by Theorcm lIII.9[ allow one to calculate the 
Euler characteristic of M. v . 

Such results have been reported in the literature for sev- 
eral model systems. As an illustrating example, we re- 
view results for the critical points and for the Euler char- 
acteristic of the mean-field /c-trigonometric model. 
This model is characterized by the potential 



Vk(q) 



A 



TV 

*1, ■■■,ik- 



[1 - cos (q i± 



Qi>)], (3-4) 



where A > is some coupling constant and N is the 
number of degrees of freedom of the system. The position 
coordinates qi € [0, 2ir) (i = 1, . . . , N) are angular vari- 
ables, so that the configuration space has the shape of an 
iV-dimensional torus. The potential describes a fc-body 
interaction where k £ N, and the interaction is of mean- 
field type, i. e., each degree of freedom interacts with each 
other at equal strength. For this model, a number of 
thcrmodynamical as well as topo l ogical quanti ties have 
been computed bv lAngelani et al\ (|2003l l2005af ). the lat- 
ter ones by making use of Morse theory. 11 Among those 
results we mention the following: 

(i) In the thermodynamic limit N — > oo, the mean- 
field fc-trigonometric model does not show a phase 
transition for k = 1, whereas it has a phase tran- 
sition for k ^ 2. The transition is continuous for 
k = 2 and discontinuous for all k ^ 3. The poten- 
tial energy at which the phase transition occurs is 
vt = A, where the index t is for transition. 

(ii) For any finite N, the ith component of the critical 
points of the potential \% is given by 



k- 1 



= 7r < m, - 



N 



N 



(mod 2tt), (3.5) 

where £ € Z, m = (mi, . . . , toat) G {0, 1} , and 
O denotes the Heaviside step function. For our 
purposes it is important to notice that the num- 
ber of critical points increases unboundedly with 
N. From the critical points and their indices, the 
Euler characteristic x (-M v ) can be computed. 



11 Vfe as given in Eq. 1 13. 41 1 is not a Morse function, b ut, in the spirit 
of Eq (TJ2J, it can be perturbed into one; see lAngelani et al\ 
(20053). 
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(iii) In the thermodynamic limit, the critical values 

1 



,,m,e 



N 



(3.6) 



become dense on the interval [0, 2 A]. This leads to 
a continuously varying (on [0,2 A]) limiting distri- 
bution 



a{v) = Jim —]n.\x(M v )\ 

N— *oo iv 

= —n(v) lnn(v) — [1 — n(v)] In [1 — n(v)] 



(3.7) 



of the modulus of the Eulcr characteristic x {-Mv) 
of M. v , where 



n(v) 



1 






1/fc" 


2 


l-=gn(l-£) 







(3.8) 



[see lAngelani e£ aLl (2005a) for the derivation of 
this result]. It is a remarkable observation that a, 
which is a purely topological quantity, already sig- 
nals the absence or presence of a phase transition 
(see Fig. Q] for a plot of the graph of a): In the case 
of k = 1 where the system does not show a phase 
transition, a is a smooth function. For k J? 2 there 
is a nonanalytic point of cr at i> = A, which is pre- 
cisely the value v t of the potential energy at which 
the phase transition occurs. 




FIG. 1 (Color online) Logarithmic modulus a of the Euler 
characteristic of M v as a function of the potential energy v for 
the mean-field fc-trigonometric model in the thermodynamic 
limit [see Eqs. $£7} and flU}]. 

In particular this last observation — together with similar 
findings for other models — may serve as a motivation to 
study the relation between thermodynamic phase transi- 
tions and the topology of the configuration space subsets 
M v . 



published to date ( Franzosi et all l2000f) . The key idea 
of this work is sketched in the following, and the results 
of an application to the </? 4 model on the square lattice 
are summarized. As an outlook, a different numerical 
approach is sketched for which an application is currently 
under way. 



1. Euler characteristic via Gauss-Bonnet theorem 

A remarkable theorem, found independently by Gauss 
and Bonnet, connects geometrical and topological quan- 
tities of manifolds. 

Theorem III. 10. Let M be a compact, orientable Rie- 
mannian manifold of even dimension n without bound- 
ary. Then its Euler characteristic x(M) is given by 



X(M) 



f daK(a), 



(3.9) 



where | S>™ | denotes the volume of an n-dimensional unit 
sphere and K is the Gauss-Kronccker curvature. 

The Gauss-Kronecker curvature is a measure of how 
much the normal vector of the manifold at a point a 
changes upon infinitesimal variation of tr. 12 This the- 
orem is remarkable, and useful for our purposes, since 
it expresses global topological properties of a manifold 
in terms of local geometrical ones. The locality of the 
Gauss-Kronecker curvature K allows to estimate (apart 
from a prefactor) the integral in (|3.9p by implementing 
some dynamics on M and probing K in the course of the 
dynamical evolution. 

This method has b een proposed and implemented by 
iFranzosi et al\ ( 2000l ) for the numerical computation of 
the Euler characteristic of the configuration space sub- 
sets S„ of the (p 4 model in two dimensions on the square 
lattice with nearest-neighbor interactions. This model is 
characterized by the potential 13 



ynn 



(<?) = 



JV 

E 

i=l 



1 2 



1 , 
4* 



(3.10) 



where J £ 1R is a coupling constant and the angular 
brackets denote a summation over all pairs of near- 
est neighbors on the square lattice. The first term of 
the potential is an on-site potential with the shape of a 
double well, and the second term describes a pair interac- 
tion between each degree of freedom ^ e 1R and its four 
nearest neighbors on the lattice. The superscript of V£ n 



C. Numerical computation of topological quantities 

Only one numerical study of the configuration space 
topology and its relation to phase transitions has been 



12 For a precise definition see an y introductory text on differential 
geometry, e.g., Chapter 3F of lKuhnell [120021) . 

13 This definitio n differs slightly from the one used by 
Fr anzos i et al but both models can be mapped onto each 
other by a suitable rescaling of J, V, and q. 
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serves to distinguish the ip A model with nearest-neighbor 
interactions from a similar one with long-range interac- 
tions introduced in Sec. IVI.AI For positive values J > 
of the coupling constant, a parallel orientation of neigh- 
boring degrees of freedom is energetically favorable, and 
the model is known to show a phase transition from a 
ferromagnetic phase at low temperatures to a paramag- 
netic phase at high temperatures in the thermodynamic 
limit. 

Implementing a Monte Carlo dynamics on £„, 
iFranzosi et all ( 2000f ) computed numerical estimators for 
the relative variation of the Euler characteristic x(£„) of 
the nearest-neighbor tp model (see their publication for 
more details on the implementation). Their result for a 
system of 7 x 7 lattice sites is reproduced in Fig. [5J The 



x 



o 




FIG. 2 Logarithmic modulus a of the Euler characteristic of 
as a function of the potential energy per degree of freedom 
for the nearest-neighbor ip 4 model on a 7 x 7 square lattice. 
The dotted vertical line marks the potential energy at which a 
phase transition takes place in the thermodyn amic limit. The 
solid line serves as a guide to the eye. From IFranzosi et all 
(|2000h . 

data suggest that a kink is present in x(^v) at a value of 
the potential energy per particle v which is very close to 
the value of v at which a phase transition takes place in 
the thermodynamic limit. 



2. Critical points of the potential 



Morse theory, as introduced in Sec. IIII.A| suggests an- 
other, straightforward way for the numerical estimation 
of topological quantities: Determining numerically the 
critical points and their indices of a potential v, one 
can calculate the Euler characteristic of the configura- 
tion space subsets M v by means of Theorcm lIII.91 Since 
the number of critical points in generic models is assumed 
to grow exponentially with the number N of degrees of 
freedom, this task is computationally very hard, and for 
N of order 10 2 or larger one can expect to find at best a 



small fraction of the existing critical points. 

Numerical methods for the computation of critical 
points of high-dimensional functions have been employed 
and tested in the study of the properties of potential en- 
ergy landscapes of g l assy systems, clusters, and others 
(jGrigera et all , 120021: IWalesi . |2004|) . An application of 



these methods to the study of configuration space topol- 
ogy and phase transitions is currently under way. 



IV. NONANALYTICITIES IN FINITE SYSTEMS 

Before turning to the study of the relation between 
configuration space topology and phase transitions, i.e., 
the behavior of thermodynamic functions in the thermo- 
dynamic limit, we discuss the somewhat simpler case of 
finite systems. The central question addressed in this 
section is: 

What are the differentiability properties of the thermo- 
dynamic functions sjv and /at for finite numbers N of 
degrees of freedom? 

For the configurational canonical free energy the an- 
swer is well known: /jv is a smooth function on its entire 
domain for all finite values of N ( Griffiths! Il972t ). For 
technical reasons, the microcanonical ensemble has been 
used only sporadically for explicit calculations in statis- 
tical physics, and little attention has been paid to the 
differentiability properties of the configurational micro- 
canonical entropy. In this section we investigate the rela- 
tion between nonanalytic points of sjq and critical points 
of V (which, in turn, are related to topology changes of 
the configuration space subsets M v ): 

What is the effect of a critical point q c of the poten- 
tial V on the differentiability of sn at the corresponding 
critical value v c — V (q c ) /N? 

Recently, exact model calculations have been per- 
formed which demonstrate that the microcanonical en- 
tropy sjv or the configurational microcanonical en- 
tropy s 7v of classical statistical ph ysics models can have 
nona n alytic points for fin it e N (ICasetti and Kastner . 
20061: iDunkel and Hilbertl. 120061; iHilbert and Dunkel . 



20061: IKastner and Schnet d . l2006f) . 
the case of smooth poten tials V 



as 



At least 
considered 



m 

by 



IKastner and Schnet d ( 20061 ) and ICasetti and Kastnen 
(120061) where the mean-field spherical model was stud- 
ied, nonanalyticities of sn show up precisely at the crit- 
ical values v c of V. 

The occurrence of nonanalyticities in s^r — though a 
surprise even to many researchers working in the field — 
can be anticipated from the discussion of a simple one- 
dimensional example. Considering a ip potential on the 
real line, 



V{q) = -\q 2 



(4.1) 



the calculation of the configurational microcanonical en- 
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tropy yields 



Si(v) = In 



0(-«) 



(4.2) 



This function has a nonanalytic point at argument zero, 
which is identical to the critical value of V corresponding 
to the critical point q = 0. From this one-dimensional 
example, it appears plausible that nonanalyticities show 
up also at critical values of higher dimensional potentials. 
This reasoning is con firmed by the following theorem by 
iKastner et ^T(|2007f) . 



Theorem IV. 1. Let V : G — ► R be a Morse function 
with a single critical point q c of index i in an open re- 
gion G C R^. Without loss of generality, we assume 
V(q c ) = 0. Then there exists a polynomial P of degree 
less than N/2 such that at v = the density of states 
51 at = cxp(A^sjv) can be written in the form 



n N (v) = p( v ) 



h N ,i( v ) 



x/ldet [i}y(<7c 



+ o{v N/2 ^) (4.3) 



for any e > 0. Here S)v is the Hessian of V, o denotes 
Landau's littlc-o symbol for asymptotic negligibility, and 

r (_ 1 )(w-o/2(_ u) (w-2)/2 0( _ v) for ^ ■ odd; 

X i (-l)'/2 v (N-2)/2 for • evenj 

j^(_l)('+l)/2 w (A r -2)/2 7r -l ln | u | for jy • odd; 

(4.4) 

is a universal function. 

Proof. For a proof of this theorem, the density of states is 
calculated separately below and above the critical value 
v = 0. By complex continuation it is possible to subtract 
both contributions and to evaluate the leading order of 
the difference. A detailed proof will be published else- 
w here. A r e lated , but weaker result has been announced 
bv lSpinellil (ll999f ). □ 



We see from this theorem that, for all finite system 
sizes, a critical point of V produces a nonanalyticity in 
the entropy sjv at the corresponding critical value. The 
theorem refers to a single critical point, but this condi- 
tion can be readily released by considering the union of 
regions around several critical points. 

Remarkably, the type of nonanalyticity described by 
liNA in Eq. (|4.4jl does not depend on the precise value 
of the index i of the critical point, but only on whether 
iV" and i are odd or even. One can verify that, indepen- 
dently of the three cases in Eq. (I4.4|) . hjy^, and therefore 
also the entropy sn, is precisely [(N — 3)/2j times con- 
tinuously diffcrcntiable. This result is in agreement with 



the nonanalytic behavior of the exact solution for the 
entropy of the mean-field sp herical model as reported by 
IKastner and Schnetzl ( 2006| ). In other words, sn becomes 
"smoother" with increasing number A" of degrees of free- 
dom, and already for moderate N it supposedly will be 
impossible to observe such a finite-system nonanalyticity 
from noisy experimental or numerical data. 



V. PHASE TRANSITIONS AND CONFIGURATION 
SPACE TOPOLOGY 

We have seen in the previous section that, for finite sys- 
tems, critical values of the potential V (and hence topol- 
ogy changes of the configuration space subsets M. v ) are 
directly related to the occurrence of nonanalytic points 
of the configurational microcanonical entropy sjv- F° r 
infinite systems, such a straightforward correspondence 
cannot be expected since, according to Theorem IIV.11 
the order of differentiability of sn in the presence of 
a finite number of critical points of V diverges when 
A" — > oo. Furthermore, from the results on the mean- 
field fc-trigonometric model reviewed in Sec. IIII.B1 it is 
obvious that a relation between topology changes and 
phase transitions for the infinite system case has to be 
subtler: The number of critical values v c of the potential 
energy where topology changes of Ai v occur increases 
unboundedly with the number N of degrees of freedom 
of the system and becomes dense on the interval [0, 2 A] 
in the limit A^ — > oo. Since a phase transition in this 
model occurs only at one particular value of the potential 
energy v = Vt, it is clear that there cannot be a one-to- 
one relation between phase transitions in infinite systems 
and topology changes. Nonetheless, the existence of some 
sort of relation is suggested by the singular behavior of 
the topological invariant a at v t as plotted in Fig. [TJ 
Similar results, i. e., a continuum of topology changes of 
the M. v and a nonanalyticity in a quantity characteriz- 
ing the topology in the thermo dynamic limit , were found 
for the moan- field XY model (ICasetti et a?J . ll999ll2002t 
2003). Numerical studies also indicate that the same fea- 
tures are present in the ip 4 mo del on a square l attice with 
nearest-neighbor interactions (jFranzosi et al 1 120001) . 



A. Conjectures 

From the results of such model calculations, a general, 
but somewhat unspecified "relation" between topology 
ch anges of the ./Vi^ and phase transitions was conjectured 
bv lCaiani et al\ (|l997t h Later on, this conjecture became 
known as the topological hypothesis, and its formulation 
gave leeway for different interpretations of its content. 

The results on the mean-field fc-trigonomctric model, 
the mean-field XY model, and the square lattice p 4 mod- 
el cited above might suggest the following formulation. 

Conjecture V.l. If, in the thermodynamic limit, the 
logarithmic density of the Euler characteristic x of 
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has a nonanalytic point at v = v t , a phase transition 
takes place at the potential energy vt ■ 

A weaker version of such a hypothesis might be formu- 
lated as follows. 

Conjecture V.2. A topology change within the family 
{■M v } v pji at v = vt is a necessary condition for a phase 
transition to take place at v t . 

We will see shortly that neither of these conjectures 
holds true for arbitrary statistical systems, and coun- 
terexamples disproving their general validity will be given 
in Sec. I VII Nonetheless the conjectures are of interest as 
they prove to be correct at least for certain (large and rel- 
evant) classes of systems. In fact, a theorem presented 
in the next section asserts that Conjecture IV. 21 is correct 
for a class of models with short-range interactions. 



B. Franzosi-Pettini theorem 

In 2004, Franzosi and Pettini announced the proof of 
a theorem asserting that C onjecture IV. 21 is corre c t for a 
class of short-range models ( Franzosi and Pettinil . l2004h . 
In the pr esent section we sketc h thi s result, and refer the 
reader to lFranzosi et al\ ( 2007| ) and lFranzosi and Pettinil 
(I2007D 1'or details. We start with the definitions necessary 
to specify the class of models for which the theorem holds. 

As throughout the paper, we consider systems of N 
degrees of freedom described by a Hamiltonian of stan- 
dard form (|2.2[) . We assume the potential V, defined on 
the continuous configuration space Tn, to be subject to a 
number of conditions, whose definitions are given in the 
following. 



Definition V.3. The potential V is of standard form 

N N 
1=1 — \ 

if it consists of an on-site potential <f> and a pair potential 
tjj, the latter depending only on the Euclidean distance 
||-|| of the degrees of freedom. For a lattice system the 
coefficients Cy determine the coupling strength between 
the degrees of freedom on different lattice points. For a 
fluid (nonlattice) system we typically have = 1 — S^j, 
where 5ij is the Kronecker symbol. 

Definition V.4. The potential V is of short range if 

on a lattice: the coefficients cy are nonzero only for i, j 
from a finite neighborhood on the lattice; 

in a fluid: for large x, \ip (x)\ decreases faster than x~ d , 
where d is the spatial dimension of the system. 

Definition V.5. The potential V is stable if, for any 
JV 6 N and all q 6 Tn, there exists a constant B 
such that 



Definition V.6. The potential V is said to be confining 
if M. v as defined in Eq. (|2.4j) is a compact set for all finite 
values of v. 

Physically speaking, a particle in a confining poten- 
tial cannot escape to infinity at finite energy. Simple 
one-dimensional examples of confining and nonconfining 
potentials are plotted in Fig. [3J 




FIG. 3 (Color online) One-dimensional illustrations of confin 
ing and non-confining potentials. Left: The parabola V(q) = 

oo. Right: The Morse po 



q 2 is confining, as lim V(q) 



tential V(q) = (l 



q— *±oo 

- e~ q ) 2 is nonconfining, as lim V(q) < oo. 

q — »oo 



It is well known that for potentials of standard form 
which arc stable, confining, and of short range, the ther- 
modynamic limit of thermodynamic functions like the 
configurational microcanonical entropy (|2.8p or the con- 
figurational canonical free energy (|2. 10[) exists, i.e., 



(v) = lim sjv(v) < oo 



and 



A' 



/„,(/?) = lim f N (f3) < oo 

N—KX> 



(5.3) 



(5.4) 



V(q) > —NB. 



(5.2) 



[see iRuelld ( 1969f ) for a proof of these and related results 
for even larger classes of systems]. For this class of sys- 
tems, Franzosi and Pettini proved a theorem which can 
be phrased as follows. 

Theorem V.7. Let V be a potential of standard form 
which is smooth, stable, confining, bounded below, and 
of short range. If there exists an interval [ui,i>2] such 
that, for any N larger than some constant No, the 
{■Mv} ve i vi „ i do not change topology, then in the ther- 
modynamic limit the canonical configurational free en- 
er gy Jn is at least two times continuously diffcrentiable 
in the interval ((3 (vi) , (v 2 )) , where (3 (vi) and /3 (1)2) 
are the values of the inverse temperature corresponding 
to the potential energies v\ and V2, respectively. 

Proof. The proof makes use of the fact that, in an in- 
terval free of critical values of V/N, the configurational 
entropy sat is smooth for all finite N (which is a result 
known from Morse theory). Then it is shown that at least 
the first four derivatives of sn are uniformly bounded 
above in N. This implies that Sqq is three times con- 
tinuously diffcrentiable, and therefore /oo is two times 
continuously diffcrentiable in the corresponding interval 
of inverse temperatures . A detailed proof can be found 
in lFranzosi et al. I (|2007f ). □ 
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Although this theorem does not exclude the possibil- 
ity of a discontinuity in some higher derivative of foo , one 
would assume that the extension of the theorem to higher 
derivatives, though laborious, is straightforward. Such a 
generalization to arbitrary derivatives of would assert 
the correctness of Conjecture IV.2I for the class of poten- 
tials covered by the theorem. 

Theorem IV.7I rigorously establishes a connection be- 
tween phase transitions and topology changes of the M. v 
for a certain class of systems, which from a conceptual 
point of view is a remarkable result. For potentials like 
that of a /3-Fermi-Pasta-Ulam chain, having no further 
critical points apart from their overall minimum, Theo- 
rem IV. 71 allows one to rigorously exclude the occurrence 
of a phase transition. Possibly even more important, for 
systems like the ip 4 model in Eq. (|3.10p one may deduce 
an upper bound on the transition potential energy vt 
from the presence of an upper bound on the critical val- 
ues of V. 

For the typical models of interest, the critical values 
v c of the potential become dense in the thermodynamic 
limit, in some cases on an interval, but often on the en- 
tire codomain of V/N (like in the case of the mean-field 
fc-trigonomctric model discussed in Sec. IIII.B[) . In that 
case, no interval [t>i,i>2], free of critical points in the 
sense specified in Theorem IV.71 exists, and no conclu- 
sions can be drawn from the theorem. An attempt to 
derive a stronger result which sheds light on the relation 
between phase transitions and topology changes even in 
the presence of (possibly very many) cri tical points has 
been made by Franzosi and Pettinil ( 20071 ). and in a mod- 
ified form bv lKastner et al. I (l2007h . and the hope is to un- 
derstand from these studies how topology might "act" to 
produce a nonanalyticity in some thermodynamic func- 
tion in the thermodynamic limit. 



C. Models not covered by Theorem IV. 71 



tures IV. II and IV. 21 and we look at these cases in more 
detail in Sec. EH 



The above observations, together with Theorem IV. 71 
can be interpreted in the following way: A change of the 
topology of the M v under variation of the parameter v 
is one possible "mechanism" which can lead to a non- 
analyticity in a thermodynamic function. However, the 
fact that some of the models listed in Table U are not in 
accordance with Conjectures IV. II and IV.2I indicates that 
a topology change is not the only such mechanism, and 
we get to know a second nonanalyticity generating mech- 
anism in Sec. IVI.Cl Theorem IV. 71 then suggests that, for 
the — from a physical point of view very important — class 
of systems fulfilling its assumptions, a topology change is 
the only mechanism available to cause a nonanalyticity. 



VI. LIMITATIONS OF THE RELATION BETWEEN 
PHASE TRANSITIONS AND CONFIGURATION SPACE 
TOPOLOGY 



We have seen in Sec. IV. Bl that, for a certain class of 
systems, a topology change in the configuration space 
subsets M v is necessary for a phase transition to take 
place at the corresponding energy or temperature. For 
some models however, as indicated in Table HI a phase 
transition is not necessarily accompanied by a topology 
change. These results immediately suggest the following 
question: 

Which of the restrictions on the class of systems for 
which Theorem \ V. 7| holds are mere technicalities which 
could be relaxed by more refined methods of proof, and 
which ones are really essential? 

In the present section we discuss two types of sys- 
tems which disprove the general validity of Conjec- 
tures IV. II and IV. 2\ thereby showing that at least the 
short-rangedness and the confining property imposed on 
the potential V in Theorem IV. 71 cannot be relaxed. 



Theorem IV. 71 establishes a necessary relation between 
phase transitions and configuration space topology for 
a certain class of systems with short-range interactions 
and confining potentials. Of course, this result does not 
exclude the possibility that a phase transition in a long- 
range system or in a system with a nonconfining poten- 
tial may nonetheless be related to a topology change in 
configuration space. In fact, all models for which ana- 
lytic calculations of the topology of the A4 V have been 
reported so far do not meet the assumptions of Theo- 
rem IV.71 being either long range, nonconfining, nons- 
mooth, or even some combination of these properties. 
Results from model calculations reported in the litera- 
ture are summarized in Table U and we notice that for 
several of these models which do not comply with the 
assumptions of Theorem IV.71 a relation between phase 
transitions and topology is nonetheless observed, even in 
the sense of the stronger Conjecture IV. 11 Some other 
models, however, disprove the general validity of Conjec- 



A. Long-range interactions 

Systems with long-range interactions are often ne- 
glected in standard treatises on statistical mechanics. In 
many cases they are even outside the scope of traditional 
thermodynamics because, in contrast to the case of short- 

the thermodynamic 
(1531) and 



range interactions (Definition IV.4[) 
limit of thermodynamic functions as in Eqs. 
(|5.4p does not necessarily exist. Furthermore, the stabil- 
ity of thermodynamic functions, manifest i n their convex- 
ity pro perties, is no longer guaranteed [see lDauxois et al\ 
(2002) for a review of the dynamics and thermodynamics 
of systems with long-range interactions] . This limitation 
of traditional thermodynamics is remarkable in regard 
to the importance of long-range interactions in physics, 
most notably in gravitation and in electrodynamics (at 
least in the absence of screening effects) . In this section, 
we give an example of a long-range system whose phase 
transition is not accompanied by a topology change in 
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TABLE I (Color online) Results from calculations of the topology of M v for some models (for definitions of the models, see 
the references indicated). The abbreviation m.-f. is for mean-field interactions and n. n. is for nearest-neighbor interactions. 
The columns list (from left to right) whether the respective model has a phase transition, whether the results on the topology 
of M v were obtained analytically or numerically, whether the potential of the system is short range, smooth, and confining, 
and whether the results are in accordance with Conjectures IV. 1 1 and IV. 2l 



model 


phase 
trans. 


anal./ 
num. 


short- 
range 


smooth 


conf. 


Conj. 

lv.il 


Conj. 
|V.2| 


references 


m.-t. A i 


yes 


anal. 


no 


yes 


yes 


yes 


yes 


(|Casetti et all 1999. 2002. 2003|) 


m.-f. /c-trigon. 


yes 


anal. 


no 


yes 


yes 


yes 


yes 


(jAngelam et at. 2003, 2005a) 


m.-f. spherical 


yes 


anal. 


no 


yes 


yes 


yes 


yes 


(Kastncr, 2006a; Kastncr and Schnctz, 2006; 


iRibeiro Teixeira and Stariolo, 2004) 


m.-f. p 4 


yes 


anal. 


no 


yes 


yes 


no 


no 


(Andronico et all 12004: Baroni. 120021: 
Garanin et all |2004|; lHahn and Kastnerl, 
I2005T) 


n. n. spherical 


yes 


anal. 


no" 


yes 


yes 


no 


yes 


(|Risau-Gusman et all 120051) 


Peyrard-Bishop 


yes 


anal. 


yes 


yes 


no 


no 


no 


(|AnEelani et all l2005bl; iGrinza and Mossal, 
12004) 


Burkhardt 


yes 


anal. 


yes 


no 


no 


no 


no 


(Angclani et al, 2005b; Kastncr. 2004) 


2d n. n. ip 4 


yes 


num. 


yes 


yes 


yes 


(yes) b 


yes 


(Franzosi et al., 20(33) 


ldn.n. XY 


no 


anal. 


yes 


yes 


yes 






(Casctti et al, 2003) 



a Although the pair potential of the nearest-neighbor spherical model appears to be short range, an effective long rangedness is induced 

by the spherical constraint. 
b As far as discernable from numerical data. 



configuration space. In Sec. IVI.CI we argue that it is 
precisely due to the above-mentioned convexity proper- 
ties of thermodynamic functions that phase transitions in 
long-range systems are not necessarily related to topol- 
ogy changes. 

The example which we discuss is the mean-field ip 4 
model, defined by the potential 



N 

E 

i=l 



1 9 1 4 

'2?i + 7^ 




(6.1) 



with coupling constant J > 0. The first term of the po- 
tential is an on-site potential with the shape of a double 
well. The second term of V^ nf describes a pair interaction 
of mean-field type, i.e., each degree of freedom q t e R 
is coupled to each other at equal strength, which is an 
extreme case of long-range interactions. The choice of an 
TV-dependent coupling strength —J/(2N), though ques- 
tionable from a physical point of view, guarantees the 
existence of the thermodynamic functions in the limit 
N -> co. 

This model is exactly solvable in the sense that, 
in the thermodynamic limit, thermodynamic func- 
tions like the configurational microcanonical en- 
tropy (jCampa and Ruffol . l2006l ; lHahn and Kastnerl 120051 
120061) or the configurational canonical free 



2UUta) or the conhgurational canonic al tree energy 
(iDauxois eiaH l2003t [Ovchinnikov and Onisclwkl Il990l ) 
can be expressed in terms of a maximization and a single 
integration. A continuous phase transition is found to 



take place in the mean-field (/s 4 model, and from the con- 
figurational microcanonical entropy an implicit expres- 
sion for the potential e nergy i> t at which the tra nsition 
occurs was derived by lHahn and Kastnerl ( 2006f ). For 
large values of the coupling constant J, an expansion of 
this implicit expression yields 



v t (J) = aJ 2 -(2a-\)J + 0{l) 1 



(6.2) 



with a = r(3/4)/r(l/4). For our purposes it is impor- 
tant to note that Vt increases unboundedly with J, i.e., 
the transition energy can be made arbitrarily large by 
increasing the coupling constant J. 

A study of the topology of the configuration space sub- 
sets M. v of the mean-f ield ip 4 model has first been re- 
ported by [Baron 1 (12002^ and later, indepen dently and by 
different methods, by Garanin et al\ (|2004j ). Looking for 
critical points q c satisfying dV^ nf (q c ) = 0, one can show 
that the corresponding critical values 



N 



(6.3) 



are nonpositive for all c ritical points q c of the potential 
y mf [see Appendix E of lBaronil (|2002t ) for a proof] 



Confronting the results (|6 . 2[) and (|6.3[) . i.e., the un- 
bounded growth of the potential energy vt of the phase 
transition with the boundedness from above of the crit- 
ical values v c of VJ llf , one immediately arrives at the 
conclusion that, for the mean-field ip 4 model, the phase 
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transition is not accompanied by a topology change in 
configuration space in general (i. e., not for arbitrary cou- 
pling constants J). Hence neither conjccture lV.2l nor con- 
jecture IV. II on the relation between configuration space 
topology and phase transitions hold true for this model. 
The mean-field Lp A model fulfills, apart from the short 
rangedness, all requirements of Theorem IV. 71 Hence we 
conclude that the assumption in Theorem IV. 71 of the po- 
tential being short ranged cannot be relaxed in general. 



B. Nonconfining potentials 



Particles in a confining potential (Definition IV.6[) are 
restricted to a bounded subset of the configuration space 
for any finite value of the energy. Modeling gases or flu- 

the 



ids by potentials of standard form (Definition IV.3|1 
pair interaction between particles is typically assumed to 
be nonconfining, but the addition of an on-site potential 
modeling a container renders the overall potential confin- 
ing. Nonconfining potentials are of physical interest for 
modeling fluctuations of interfaces by means of so-called 
solid-on-solid models (jAbrahaml . Il986h . The occurrence 
of a localization-delocalization transition in these mod- 
els is a consequence of the nonconfining property of the 
potential. In this section, we discuss a class of solid- 
on-solid models with nonconfining potentials, for which 
we find phase transitions not accompanied by topology 
changes in configuration space. 

The solid-on-solid models we consider are one-dimen- 
sional lattice models characterized by a potential of the 
form 



N-x 



N 



V sos {q) = J ^2 



Qi+i 



(6.4) 



with some coupling constant J > 0. For the modeling 
of physical systems, the parameter n typically has values 
n = 1 or 2. The position coordinates qi can take on val- 
ues from the positive half-line, so that the configuration 
space of the model is Tn = (R + ) JV . The pair interaction, 
in contrast to the example discussed in Sec. IVI.A| is of 
short range, being restricted to nearest neighbors on the 
lattice. The on-site potential U is a real-valued function 
with a single local minimum somewhere on its domain 
R + , and it approaches a finite value in the limit of large 
arguments, 



lim U(x) < oo 



(6.5) 



(like the nonconfining one-dimensional potential in the 
right plot of Fig. [3|) . As a consequence of Eq. (|6.5p , we 
have that the potential V^os m Eq. (|6.4[) is also noncon- 
fining. 

By means of a transf er operator technique 
( Kramers and Wannierl . Il94lt ). the thermodynamic 
limit value of the configurational canonical free energy 
/oo of the solid-on-solid model (|6.4[) can be written as 

-AM/3) = max [In A 4 (/?)], (6.6) 



where are the eigenvalues of the (well-behaved) solu- 
tions of the eigenvalue problem 



/>oo 

/ dy cxp {-\pu{x) -jp\x- y\ n - \pU{y)} $ (y) 
Jo 

= \&(x). (6.7) 



Equation (|6.7p can be transformed into a Schrodinger- 
type eigenvalue problem. This m apping is exact for 
the case n = 1 ( Burkhardil Il98lh . whereas it involves 
a grad ient-expansion app r oxima tion for the case n = 
2 [see iTheodorakopoulosI (|2003l ) for details]. For the 
Schrodinger-type equation thus obtained, one can argue 
that, depending on the values of J and /3, bound state 
solutions may or may not exist, and it is this changeover 
which corresponds to the occurrence of a phase transition 
in the statistical mechanical model. For simple box-type 
potentials U, the eigenva l ue problem can be solved ana- 
lytically (|Angelani et all . l2005bl ; iBurkhardtl . Il98lh . and 
for these cases the model is found to undergo a phase 
transition at some inverse temperature Pt(J) which is a 
continuously varying, nontrivial function of the coupling 
constant J. Similarly, for arbitrary on-site potentials U, 
a nontrivial dependence on J of the transition (inverse) 
temperature (3 t as well as of the corresponding transi- 
tion potential energy v t is expected (and corroborated 
by numerics). 



From the point of view of configuration space topology, 
the class of solid-on-solid models (|6.4|) is particularly sim- 
ple. A complete charact erization of the topology of the 
M v has been outlined bv lKastnerl (|2004h , and details are 
given in the Appendix. 14 For any on-site potential U 
which is a monotonous function on the intervals (0, aWn) 
and (x m i n , oo) that are left and right to the location x m i n 
of the unique minimum of U, one finds that two topology 
changes take place within the family {M v } ve j i of con- 
figuration space subsets. One is located at the ground 
state energy v± — U (x mm ), while a second one appears at 
V2 = lim U(x). Both, v\ and vi are independent of the 

value of the coupling constant J. Comparing this finding 
with the nontrivial dependence of the transition poten- 
tial energy u t as argued above, we observe that phase 
transitions in the class of models defined by Eq. (|6.4p arc 
not related to topology changes in configuration space. 



Setting n = 2 in Eq. (|6.4|) , the potential V S os is smooth 
and the model fulfills, apart from the confining property, 
all requirements of Theorem IV.7I We therefore conclude 
that the assumption in this theorem of the potential be- 
ing confining cannot be relaxed in general. 



14 This r esult includes, as a special case, those of[Grinza an d Mossal 
(2004) for the Peyrard-Bishop model, obtained by a different 
method. 
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C. Nonanalyticities from maximization 



The examples presented in Sees. [VO] and [VLB] dem- 
onstrate that, at least for systems with long-range in- 
teractions and for systems with nonconfining potentials, 
a phase transition is not necessarily accompanied by a 
topology change in configuration space. In these cases 
one might suspect, instead of a topology change, a dif- 
ferent kind of mechanism behind the occurrence of a 
phase transition. In this section we argue that for both, 
the cases of long-range systems and nonconfining po- 
tentials discussed in Sees. IVI.AI and I VLB! nonanalyt- 
icities in thermodynamic functions are generated from 
smoot h functions by means of a maximization mecha- 
nism (|Kastneij . l2006bT ). 



1. Mean-field </? 4 model revisited 

The long rangedness of the pair potential of a system 
has, as mentioned in Sec. IVI.AI remarkable consequences 
on the convexity properties of the thermodynamic func- 
tions. Notably, the (configurational) microcanonical en- 
tropy of a system with long- range interactions, in con- 
trast to the short-range case dGallavotti . 19991; Lanford . 
1973), is not necessar ily concave ( Dauxois et all 2002 : 
Touchette et a/.l . l2004h . 



Definition VI. 1. (i) A set A C TR™ is called a convex 
set if ax + (1 — a)y <G A for all x,y E A, a e [0, 1]. 

(ii) A function / : A — > JR defined on a convex set A is 
called a convex function if 

f{ax + (1 - a)y) ^ af(x) + (1 - a)f(y). (6.8) 

(iii) If — / is convex, / is called a concave function. 

In the case of long-range interactions, a nonana- 
lytic point of the configurational microcanonical entropy 
s»D (v) can arise from the maximization over one variable 
of a smooth, but nonconcave entropy function Soc(v,m) 
of two variables (defined shortly), and this is precisely 
what happens for the mean-field (p 4 model. 

For the mean-field ip 4 model defined by the potential 
(|6.ip . an exact calculation of thermodynamic functions 
is possible, and large deviation techni ques are an elegant 
way to perf orm such a calculation. lHahn and Kastnerl 
( 20051 . l2006h reported exact results for two related mi- 
crocanonical thermodynamic functions of the mean-field 
(p 4 model in the thermodynamic limit: for the configu- 
rational microcanonical entropy Soo(v) as defined in Eqs. 
(|2.8p and (|5.3[) , and for the configurational microcanoni- 
cal entropy 

Soo(w,to) 

i r N 

= lim — In / dqS [V{q) - Nv] S [ 2J ft - Nm] 

N^co N J Fn j=l 

(6.9) 



as a function of two variables, namely the potential en- 
ergy v and the magnetization to. The entropy Sqq is 
found to have a nonconcave part ( Hahn and Kastnerl . 
I2005T) . Furthermore, as a consequence of a general 
result from l a rge d eviation theory [Theorem 1.4 in 
Iden Hollander! (|2000l )]. is a smooth function in both 
its variables. Applying Laplace's method for the evalu- 
ation of asymptotic integrals to Eq. (|6.9p . one can show 
that the entropy functions Soo(v) and Soo(v,m) are re- 
lated by 



(v) 



max \s. 



,(v,m)] 



(6.10) 



The maximization over m of a smooth and nonconcave 
entropy function Soo(v,m) then may lead to a nonana- 
lyticity in Soo- 

To illustrate how a nonanalyticity may emerge from a 
maximization over one variable of a smooth but noncon- 
cave function, it is instructive to consider the following 
simple example. The function g~i : TR 2 — > TR with 



gi(v,m) = v 



2vm 



(6.11) 



is, as is easily verified, a nonconcave function, and its 
graph is shown in the plot of Fig. [4] (top). Maximizing 
with respect to the second variable of g~\ , one obtains the 
function 



gi(v) = max g x (v, to) 



V a fOTW<0 ' (6.12) 
v — V for v ^ 0, 



which has a nonanalytic point at v = 0. In contrast, for 
a concave function like 



<72(i>, m) = v — v 2 



2m z 



(6.13) 



[Fig. [?] (bottom)], a maximization with respect to m 
yields 



.92 (v) = max 92 (w, m) = v — v 2 



(6.14) 



which is smooth on TR. 

Similar to the simple examples of the functions g~\ and 
gi, the nonanalytic point of the entropy Sqo of the mean- 
field y 4 -model is created from a smooth but nonconcave 
entropy Sq^i^to) by a maximization with respect to m 
(for a plot of the graph of Soo(v, to), see Fig.[S]). We inter- 
prete this maximization as the nonanalyticity-generating 
mechanism which is at the basis of the phase transition 
of the model. 

In the thermodynamic limit (if it exists), the mean- 
field approach is exact for many systems wit h long-range 
interactions [see iBouchet and Barrel (|2005h and refer- 
ences therein for a detailed discussion] . Since the mean- 
field variational problem can be written as a maximiza- 
tion of a smooth many-variable function, the same can 
be done for the entropy of these long-range systems. The 
different types of phase transitions which may result from 
such a maximization depend on the number of variables 
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FIG. 4 (Color online) Plots of the graphs of the nonconcave 
function gi(v, m) (top) and of the concave function g2(v,m) 
(bottom). 



sponding reduced entropy function of n — 1 variables, 

^ {x\ , . . . , Xi— i , Xi+i , . . . , x n ) 

= max[s(") (*!,...,<)], (6.15) 
is again smooth and concave. 

Proof. By elementary calculus or geometric considera- 
tions. □ 

For systems with short-range interactions, this re- 
sult rules out the occurrence of nonanalyticitics from 
the above described maximization mechanism: If sW is 
smooth, then, as a consequence of its concavity, s(" _1 J 
has to be smooth as well and no phase transition takes 
place. Nonetheless, in short-range models a different type 
of maximization can be the origin of a nonanalyticity, and 
this is the content of the following section. 



2. Solid-on-solid models revisited 

For the solid-on-solid models defined by potential (|6.4[) . 
the configurational canonical free energy foa is expressed 
in Eq. (|6.6| as the logarithm of the largest of the eigen- 
values Xi of the eigenvalue problem (|6.7p . 

- / 5/oo(/9) = max[lnA i (/3)], (6.16) 




FIG. 5 (Color online) Plot of the graph of the entropy func- 
tion Soo(w,m) of the mean- field (p 4 model with coupl ing con- 
stant J = 1. Adapted from lHahn and Kastnerl (|2006T ). 



and the nonanalyticities of /oo are discussed by studying 
the behavior of A,;. The Xi are expected to be smooth 
functions of f3, but they may have crossing points. 15 
If we assume this claim to be correct, the maximiza- 
tion over the index i in Eq. (|6.16j) in the presence of 
a crossing point of the largest and the second largest Xi 
is the only remaining source of a nonanalyticity in /oo- 
For short-range models like the solid-on-solid models dis- 
cussed here, equivalence of statistical ensembles holds. 
As a consequence, the entropy Soo can be obtained from 
the free energy /oo by means of a Legcndre-Fenchel trans- 
form, and the nonanalytic point in /oo gives rise to a 
nonanalytic point in Soo. 



D. Two nonanalyticity generating mechanisms 



of the entropy function and on the symmetries present. 
A complete classification of the various nonanalyticities 
which can occur for an entropy function of two var iables 
has been worked out bv lBouchet and Barrel ( 20051 ). 

In models with short-range interactions, the above- 
described mechan i sm cannot occu r due to the concavity 
(|Gallavotti Il99fl iLanfordl . Il973h of entropy functions. 
This is the content of the following theorem. 

Theorem VI. 2. Let sW(ii, . . . ,x n ) be a smooth, con- 
cave entropy function of n variables. Then the corrc- 



Two counterexamples disproving the general validity of 
Conjectures IV. II and IV.2I on the relation between phase 
transitions and configuration space topology were pre- 
sented in Sees. IVI.AI and IVI.BI For these models — 
the mean-field ip model and a solid-on-solid model — we 
have shown that nonanalyticities in the thermodynamic 



15 Related results exist in the perturbation theory of linear oper- 
ators. Modifying those to turn this claim into a rigorous result 
would be a worthwhile task. 
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functions may be viewed as arising from the maximiza- 
tion over one variable (or one discrete index) of some 
smooth function. This maximization can be interpreted 
as one possible mechanism generating a nonanalyticity, 
whereas (certain) topology changes in the configuration 
space subsets M. v are another such mechanism. In prin- 
ciple, in a given system any of these mechanism, or even 
both, may occur and trigger the occurrence of a phase 
transition. Theorem IV. 71 then asserts that, for the class 
of (short-range, non-confining, etc.) systems fulfilling its 
assumptions, topology changes are the only nonanalytic- 
ity generating mechanism at ones disposal. 

The type of nonanalyticity generating mechanism oc- 
curring has notable consequences also for physically rel- 
evant quantities: In case of a continuous phase tran- 
sition, one can show that generically 16 the critical ex- 
ponents characterizing the nonanalytic point take on 
mean-field values whenever the nonanalyticity was cre- 
ated by a maximization a s in Eq. (|6.10p [sec Appendix of 
lHahn and Kastnerl ( 20051 ) for details]. 

The identification of a second nonanalyticity generat- 
ing mechanism does by no means diminish the interest 
in the topological approach to phase transitions. The 
typical systems of interest in statistical mechanics have 
short-range interactions and confining potentials, and an 
investigation of the relation of phase transitions and con- 
figuration space topology in such systems is of great in- 
terest. However, for this class of systems analytic calcu- 
lations of topological quantities seem out of reach, which 
underlines the importance in further development of nu- 
merical techniques as in Sec. IIII.C1 



gain insights into the fundamental mechanisms which are 
at the origin of a phase transition. No final answer to 
this question has been given so far, but some of the 
preliminary results are worth mentioning. The follow- 
ing sections discuss several proposals which have been 
suggested, more or less explicitly, on the basis of the few 
model calculations which are available. 



A. Simultaneous attachment of O(N) different handles 



ICasetti et all ( 20031 ) computed critical points and in- 
dices of the configuration space subsets Ai v of the mean- 
field XY model with and without an external magnetic 
field. This model is characterized by the potential 



Vxr(g) = — [l-cos(g i -g,-)]-' l X) COS9i ' (7-1) 



2N 

where J > is a coupling constant, h € R is an external 
magnetic field, and the coordinates £ [0, 2tt) are angu- 
lar variables. The potential energy vt at which a phase 
transition occurs in this model for h = is found to 
coincide with the only critical value at which the topol- 
ogy of the A4 V involves the simultaneous attachment of 
handles 17 of O(N) different types. This observation led 
ICasetti et all (|2003l ) to the following conjecture. 



Conjecture VII. 1. A topology change of M. v at some 
V = v t which involves the simultaneous attachment of 
handles of O(N) different types entails a phase transition 
at v t . 



VII. SEARCH FOR A SUFFICIENCY CRITERION 

We turn our attention back to the class of short-range 
systems for which, according to Theorem IV. 71 a topol- 
ogy change within the family {-M v } vGR of configuration 
space subsets is necessary for a phase transition to occur 
at the corresponding energy or temperature. Although 
this theorem indicates that some sort of relation between 
phase transitions and configuration space topology ex- 
ists, it does not have much to say about the form of this 
relation. As pointed out in Sec. [Vj a topology change of 
M v is not sufficient for a phase transition to take place, 
and the obvious question to ask is: 

Under which conditions do topology changes give rise 
to a phase transition ? 

This search for a sufficiency criterion, specifying the re- 
lation between phase transitions and topology changes, 
may be considered as the big open question in the field, 
and from an answer to this question one can expect to 



However, a counterexample to this conjecture can be 
constructed by considering the mean-field XY model in 
the presence of an external magnetic field h. In this case, 
the model docs not show a phase transition, but the same 
type of topology changes as in the absence of a field h. 



B. Nonanalyticities of the Euler characteristic 



In the same paper ( Casetti et all [20031 ) . a nonanalytic 
point in the logarithmic density of the Euler character- 
istic x(-M. v ) of the mean- field XY model at v = Vt is 
reported, and the same feature is observed for the mean- 
field fc-trigonometric model (see Sec. IIII.Bl and Fig. [TJ. 
These findings suggest the following conjecture. 

Conjecture VII. 2. A nonanalyticity at v = vt in the 

(logarithmic density of the) Euler characteristic of M. v 
entails a phase transition at v t . 



Meaning "whenever none of the leading coefficients of a Taylor 
expansion of s in the vicinity of the phase transition point acci- 
dentally vanishes." 



See T heorem HTl. 71 of the present paper or Sec. 3 of iMatsumotol 
(2002) for the definition of a handle in topology and for the 
construction of manifolds by the attachment of handles. 
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This is a "sufficiency version" of Conjecture IV. 11 but 
it suffers from the same shortcoming as the — presumably 
related — Conjecture IVH.1I For nonzero external mag- 
netic field, the nonanalyticity in the Euler characteristic 
of the mean-field XY model persists, although no phase 
transition is present. 



C. Nonpurely topological sufficiency conditions 

For two different models (i. e., with and without exter- 
nal magnetic field), we have seen that in the presence of 
the same kind of topology changes the thermodynamic 
properties can differ dras tically, and this behav i or can be 
understood as follows [see lCasetti and Kastnerl ( 20071 ) for 
a related discussion] . 

The potential Vxy /N of the mean-field XY model is 
a bounded above and below function, and the value vt 
at which O(N) different handles are attached simultane- 
ously (or at which the Euler characteristic has a nonan- 
alytic point) equals its upper bound, 



1 

N 



sup Vxv{q)- 



(7.2) 



In the absence of an external magnetic field, the config- 
urational microcanonical entropy sat is a monotonously 
increasing function, and its slope has a positive lower 
bound, 



inf^M>0 

v OV 



(7.3) 



[see Fig. |6] (top) for a plot of the graph of Soo . The point 
Vt at which O(N) different handles are attached corre- 
sponds to the "end point" of sn at the upper bound- 
ary of its domain. Switching, via Eq. (|2.14[) . to the 
canonical ensemble by considering the Lcgendrc-Fcnchcl 
transform of sat, this end point is conjugate to the value 
0i = mf v [ds]\[(v)/dv} of the inverse temperature, and 
it is at this value that a phase transition occurs in the 
mean-field XY model. 

In contrast, for nonzero external field h, even if arbi- 
trarily small, the shape of s;y changes drastically. The 
graph of sat then looks as plotted in Fig. [6] (bottom), and 
its slope ds]y(v) / 'dv is not anymore bounded below. The 
value vt at which O(N) different handles are attached 
corresponds to the maximum of sjy. Canonically, how- 
ever, the temperature at v t is infinite 



T 1 - 



ds N {v) 
dv 



(7.4) 



Therefore the corresponding macrostate is thermody- 
namically not accessible and the associated "strong" 
topology change does not affect the thermodynamic be- 
havior of the system. 

Taking into account the above considerations, a suffi- 
ciency criterion of purely topological nature for the exis- 
tence of a phase transition is unlikely to exist. "Strong" 
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FIG. 6 Configurational microcanonical entropy Soa(v) of the 
mean-field XY model in the thermodynamic limit. In the 
upper plot, which is for zero external magnetic field, the slope 
of Soo has a positive lower bound. For an external magnetic 
field h = 1/2, the graph of is shown in the lower plot, 
and the slope of Soo is unbounde d above and below. Adapted 
from lCasetti and Kastnerl (120071 ). 



topology changes in the sense of conjecture IVII.ll or 
IVII.2I [attachment of handles of O(N) different types or 
a nonanalyticity of the Euler characteristic] are reason- 
able candidates for being part of a sufficiency criterion, 
but — as follows from the discussion of the mean-field XY 
model with external magnetic field — apparently have to 
be supplemented by a nontopological condition (presum- 
ably comprising some notion of measure on phase space 
or configuration space). 18 Note, however, that in gen- 
eral the influence of "measure" may be of a more subtle 
kind than the "infinite-temperature" -argument applying 
to the mean-field XY model. 



18 A sufficiency criterion for the existence of a phase transition, 
consisting of a topological part a n d a p robabilistic part, has 
been given by [Baroni an d Casettil ([2006 ) for discrete symme- 
try breaking phase transitions. The probabilistic ingredient of 
their theorem is a Peierls-type argument, and the phase transi- 
tion to which the criterion refers does in general not take place 
at the energy of the topology change considered, so — although 
makin g reference to topology in configuration space — the crite- 
rion in lBaroni and Ca sctti (2006) is not quite in the spirit of the 
topological criteria discussed before. 
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VIII. SUMMARY 



IX. EPILOG 



We have reviewed and critically discussed a topologi- 
cal approach to phase transitions, investigating the rela- 
tion of phase transitions in classical statistical mechanical 
systems and topology changes of the configuration space 
subsets M. v as defined in Eq. (|2.4[) . For the computation 
of the topology of the A4 V , the mathematical framework 
of Morse theory is particularly convenient, and we have 
summarized some of its elementary results in view of this 
application. 

In finite systems, topology changes of A4 V are inti- 
mately related to nonanalytic points of the configura- 
tional microcanonical entropy sjy. Any critical point of 
a potential V gives rise to a nonanalytic point of sn(v) 
at v = V(q c )/N. The form of the nonanalyticity is spec- 
ified in Theorem IIV.1| and the order n of the derivative 
for which d n SN/dv n becomes discontinuous increases lin- 
early with N. 

In the thermodynamic limit, the relation between non- 
analytic points of thermodynamic functions (i.e., phase 
transitions) and topology changes of M. v is more intri- 
cate. For some class of short-range systems with smooth, 
nonconfining, and bounded below potentials, a topology 
change of the subsets A4 V at v = Vt is necessary, but not 
sufficient for a phase transition to take place at v t . In con- 
trast, in systems with long-range interactions or systems 
with nonconfining potentials a phase transition need not 
be accompanied by such a topology change, as demon- 
strated by means of two counterexamples: the mean-field 
(and therefore long-ranged) ip 4 model and some class of 
solid-on-solid models with nonconfining potentials. For 
such systems, the nonanalytic point in a thermodynamic 
function can be viewed as emerging from a maximization 
over some smooth function. 

In summary, two different mechanisms which may 
cause a nonanalyticity in a thermodynamic function have 
been identified: first, certain topology changes within the 
family {-M.v} ve ]s. of configuration space subsets; and, sec- 
ond, a maximization over one variable of a smooth func- 
tion of several variables. Theorem IV. 71 then asserts that 
only the former one of these mechanism can occur in the 
class of short-range systems which are in accordance with 
the theorem's assumptions. This is a remarkable finding, 
since this class of systems contains the types of the sys- 
tems which are typically of interest in statistical physics. 

It remains an open task to precisely specify which 
topology changes entail a phase transition. Several pro- 
posals for conditions on topology changes of the A4 V , al- 
legedly sufficient to guarantee the occurrence of a phase 
transition, are discussed, but a final answer to this ques- 
tion is still lacking. One may conjecture that such a crite- 
rion will not be exclusively of topological character, but 
instead may involve some notion of measure or geometry 
as well. A solution to this problem will be a major step 
forward towards an understanding of the origin of phase 
transitions in classical statistical mechanical systems. 



From the above discussion, it should be obvious to the 
reader that the relation of phase transitions to topology 
changes in configuration space is not a settled issue, but a 
topic of current research activity. Despite its incomplete 
status, one may profit from the study of the results to 
date in various ways. 

First, the topological approach stimulates the study of 
phase transitions from a viewpoint quite different from 
the conventional one, and such a change of perspective 
may help to deepen the understanding and to inspire 
further research activity. In particular for the study of 
the relation between phase transitions and the chaoticity 
of the underlying dynamics of the system, the tqpolog - 
ical viewpoint may be beneficial ( Casetti et all l2000f h 



Again from a conceptual point of view, the topological 
approach, at least within the framework of Morse theory 
as used throughout the present paper, has remarkable 
similarities to the study of glassy systems, biomolecules, 
or clusters from t he sad d le po ints of their potential en- 
ergy landscapes (jWaled . 120041 ). and one may hope to 
profit from these similarities of methods in future inves- 
tigations. 

Apart from such general and conceptual considera- 
tions, early efforts have been seen to bring to fruit 
differential geometrical concepts, and possibly also the 
related topological c oncep ts, in physics applications. 
M azzoni and Casettil (I2006D considered curvature fluctu- 
ations of constant potential energy submanifolds of mini- 
malistic models of proteins, finding that good folders may 
be distinguished from bad folders by studying the energy 
dependence of these fluctuations. A connection between 
curvature fluctuations and t opology chang e s in c onfigu- 
ration space as described bv lCasetti et all ( 2000h might 



then lead to a topological interpretation of this result. 

Last but not least, presenting and advocating a con- 
cept applicable to classical statistical mechanical systems 
inevitably provokes the question of how an extension to 
quantum mechanical systems might look. The obvious 
idea would be to study the topology of constant en- 
ergy subsets of the underlying Hilbert space (or of some 
related projective space), but to the knowledge of the 
author no noteworthy effort along these lines has been 
made so far. What may be of related interest is a study 
of the topology of constant energy surfaces in the com- 
plex projective space o f pure qua ntum states r e porte d by 
iBrodv and Hughstonl (|200lD and iBrodv et al\ (|2007f) for 
quantum mechanical one-particle systems, and one might 
take these results as a starting point for a future investi- 
gation of phase transitions in quantum systems and their 
relation to state space topology. 
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APPENDIX A: Configuration space topology of 
solid-on-solid models 

In this Appendix, a theorem characterizing the topol- 
ogy of the configuration space subsets M. v for the class 
of solid-on-solid models (|6.4p defined and discussed in 
Sec. IVI.Bl is proven. 

Theorem A.l. Consider the potential 



N-l 



TV 



v( q ) = jJ2 h+i - r + E u (ft) ( A1 ) 

4=1 i = l 

with J > 0, n G N, and $ £ R+ for all i = 
1,..,,N. The on-site potential U : R + -> 1R is 
supposed to have a unique minimum for some argu- 
ment x mnl , and we define U mnl = U (x nim ). Let U be 
monotonously decreasing on the interval (0, a; m i n ) and 
monotonously increasing on (cc m j n ,oo). Furthermore, 
let lim U{x) = Uoo < oo, which makes V a noncon- 

x — >oo 

fining potential. Under these conditions, the subsets 
M v = (R + ) N | V (q) < vN} fulfill the equivalence 
relations 

{0 for v < U min , 

I N for U min < v < U^, (A2) 

R+ x I^" 1 for Uoo < v, 

where ~ denotes homcomorphicity, is the empty set, 
and I = [0, 1] the unit interval. 

Proof. In a first step, M. v , if not empty, is shown to be a 
star convex subset of (R+) , i. e., there exists age M v 
such that the line segment from q to any point in A4 V is 
contained in M. v . This is proven by observing that, for 
every q = ...,<?jv) 6 (R + ) Af and every A G [0, 1], the 
inequality 



V (A {q - <7min)) 

A" J\q i+ i - q t 



N-l 

E 

i=l 



N 



E U (ft ~ ^min)) 



(A3) 



s$ J|^+i - g,- 

< V - gmin) 



^ u( qi 



i) 



holds. Here <7 m i n = (x m i n , . . . ,x m i n ), and the inequality 
for the second term in Eq. (|A3[) is a consequence of the 
monotonicity properties of U. The star convexity of M. v 
implies homotopical equivalence to I N (or to an A^-ball 
B w ), but not necessarily homeomorphicity. 

In a second step, the (un)boundedness of Ai v is in- 
ve stigated. This is done, a nalogously to the treatment 
bv lGrinza and Mossal ( 2004f ). by studying the asymptotic 
behavior of V(Xq) in the limit A — > oo. As a consequence 
of the nonconfining character of U , we find 



lim V(Xq) = 



NUoo < oo if qi 
oo else. 



qj V«, j = 1, 



N, 



(A4) 



Hence for v = V/N < U^, only configurations q £ 
(R + ) JV with finite (Euclidean) norm are accessible, 
whereas configurations of arbitrarily large norm can be 
attained for v > U^. From this observation it follows 
that AA V is a bounded subset of (R + ) W for v < Uoo, 
and, together with the star convexity shown above, M v 
is found to be homeomorphic to 1^. For v > Uoo, how- 
ever, M. v is unbounded. Since configurations of arbitrar- 
ily large norm can be attained only "in a single spatial di- 
rection" , i. e., in the vicinity of the (hyper)space diagonal 
q = A(l, 1), A > 0, we conclude that M% is topologi- 
cally equivalent to the product M% ~ R+ x I^ -1 . With 
the immediate observation that M. v = for v < U m i n , 
the proof of Theorem I A. II is complete. □ 
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